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We study a frustrated 3D antiferromagnet of stacked Ji — J2 layers. The intermediate 'quantum 
spin liquid' phase, present in the 2D case, narrows with increasing interlayer coupling and vanishes 
at a triple point. Beyond this there is a direct first-order transition from Neel to columnar order. 
Possible applications to real materials are discussed. 
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I. INTRODUCTION 

The study of frustrated quantum antiferromagnets re- 
mains an active field, characterized by strong interplay 
between theory and experiment. A model which has been 
studied extensively (see [iHIBl and references therein) is 
the so called 'Ji — J2 model', a spin-1/2 system on the 
two dimensional square lattice with nearest and next- 
nearest neighbor interactions of strength Ji and J2 re- 
spectively, both being antiferromagnetic. The Neel or- 
der at T = 0, which pertains for J2 = 0, is destabilized 
by the frustrating J2 interaction and vanishes at around 
J2/J1 ~ 0.4. In the opposite limit, of large J2, a colum- 
nar (often ambiguously termed 'collinear') ordered phase 
occurs, in which successive columns (or rows) of spins 
alternate in direction. The columnar phase becomes un- 
stable at J2/J1 — 0.6. A magnetically disordered region 
thus exists in the region 0.4 ^ J2/J1 0.6, the nature of 
which remains not fully resolved. 

It is of interest to ask what happens to this intermedi- 
ate phase, and indeed to the entire phase diagram, when 
these Ji — J2 layers are coupled, perhaps weakly, in the 
third direction, forming a 3-dimensional structure. This 
is the main topic of this paper. It has already been stud- 
ied by Schmalfuss et al. [1^ using the coupled-cluster 
and rotation-invariant Green's function methods. They 
found that the disordered region becomes narrower as a 
function of J2 when J3 is increased, and vanishes entirely 
at J3/J1 ~ 0.2 — 0.3. We will address the question using 
series expansion methods in both Neel and columnar 
phases, as well as first order spin- wave theory. 

It has been argued recently that the layered materi- 
als Li2VOSi04, Li2VOGe04 are well represented by the 
spin-1/2 J1 — J2 model with J2 3> Ji, i.e. in the columnar 
phase. In ref. (l8| high-temperature expansions for the 
specific heat and magnetic susceptibility for the purely 
2-d model were used to try to constrain the values of the 
exchange parameters by fitting to the experimental data. 
At the same time an ab initio local density approxima- 



tion (LDA) calculation of the exchange parameters Ji, 
J2, J3 yielded (0.75K, 8.8K, 0.25K) for the Si material 
and (1.7K, 8. IK, 0.19K) for the Ge system. This suggests 
that the coupling in the third direction is by no means 
negligible and ought to be included in a fitting procedure. 

We mention also that a popular, though not universally 
accepted, scenario to understand the magnetic properties 
of the recently discovered superconductin g iro n pnictides 
is via a spin-1 layered Ji — J2 model [T9lT20j. While we 
do not consider these systems explicitly here, our results 
may have some relevance. 

Our model is a 3-dimensional spin-1/2 antiferromag- 
net, on a tetragonal lattice, as shown in Figure [IJa). 
Frustrating J2 interactions occur in the x-y plane but 
not in the third direction. All interactions are antiferro- 
magnetic. 

The Hamiltonian, in standard notation, is 



H 



(1) 

Ji ^ Si.Sj 



(2) 



(3) 

Ji ^ Si.Si 



(1) 



<ik> 



where the summations are over the three classes of cou- 
pling respectively. 

In the following section, section |TT1 we present our zero 
temperature calculations and results, for the ground state 
energy and sublattice magnetization (order parameter), 
for both the Neel and columnar phases. Our results 
support a phase diagram of the form shown in Figure 
[TJb), and we estimate the location of the triple point 
at J2/J1 0.54(3) and J3/J1 = 0.16(3) Our series re- 
sults are compared with the results of linear spin-wave 
theory. In section IIIII we use series methods to compute 
the 1-magnon dispersion curves, in both ordered phases, 
and again compare our results with the spin- wave predic- 
tions. The spin-wave theory follows standard lines and, 
for completeness, is outlined in the appendix. In section 
IIVI we present our conclusions and discussion. 



II. GROUND STATE BULK PROPERTIES 
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The series expansion method is based on perturbative 
calculations for a sequence of finite connected clusters, 
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(a) 

Figure 1: (a) Tetragonal unit cell of the model; (b) Conjectured 
(probably) first-order and the solid circle is thus a triple point. 

which are then combined to obtain a series for the bulk 
system. In practice it is possible to treat of order 10^ dif- 
ferent clusters and to obtain series of order 10-20 terms, 
depending on the model. The Hamiltonian is written in 
the standard form H = Hq + XV where Hq has a sim- 
ple known ground state. In the present work we use an 
Ising expansion in which Hq consists of the diagonal 
terms and the quantum fluctuations are included in the 
perturbation V. Series are obtained in powers of A, and 
extrapolated to A = 1 via standard Pade or differential 
approximant methods. The interested reader is referred 
to the book [l3| for more detail. 

A. Neel Phase 

It is well known that the unfrustrated square lattice 
(J2 = J3 = 0) has reduced Neel order in the ground 



We have added and subtracted constant terms to make 
the unperturbed energy zero. 

We have obtained series for the ground state energy Eq 
and the magnetization M to order A^. The rapid prolif- 
eration of clusters with 3 bond types (there are 320274 
with 9 or fewer sites) limits the length of the series ob- 



phase diagram of the model at T = 0. The transition lines are 



state. This is a 2-sublattice structure with spins pointing 
in opposite directions on the two sublattices. For tech- 
nical reasons it is useful to perform a spin rotation on 
the B sublattice, making the unperturbed ground state 
ferromagnetic. The Hamiltonian, for a lattice of N sites 
is then written as 



H = -^{2,h-2J2 + J3)N + Ho + XV (2) 



with 



(3) 



(4) 

I 

tainable. The data are far too extensive to present, but 
can be provided on request. The values of Eq and M can 
then be estimated, with some uncertainty, for any values 
of J2 and J3. It is convenient to set Ji = 1 to set the 
energy scale. A display of the results is deferred to the 
third subsection, where they are presented together with 



(1) 
<»j> 



(2) 



Ho = Jiy S^S]) + .h J2 i^^^k - 2 ) + ^3 > ' ( - - SfS, 



<ik> 



(3) 
<il> ^ 



, (1) -, (2) (3) 

^ = 9-^1 E ^^tst + S-S-) + - J2 ^ {sts^ + s-st) + - J3 E ^^tst + s-s^ 



<ik> 
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Figure 2: Series estimates of the ground-state energy per site at couplings (a) J3/J1 = 0.0, (b) J3/J1 — 1.0. Filled circles - 
Neel phase; fiUed triangles - columnar phase. The dashed lines are linear spin-wave predictions in each phase. 
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Figure 3: Series estimates of the staggered magnetizations in each phase at couplings (a) J3/J1 = 0.0, (b) J3/J1 — 1.0. Filled 
circles - Neel phase; filled triangles - columnar phase. The dashed curves are linear spin-wave predictions in each phase. The 
dashed vertical line in Figure |3}3 marks the estimated position of the first-order transition. 



the results for the columnar phase. 



B. Columnar Phase 

In the columnar phase the spins on alternating columns 
in a particular x-y plane will point in opposite directions. 
This satisfies all of the strong J2 bonds but leaves half 



of the Ji bonds frustrated. In adjacent planes these are 
shifted by one lattice spacing, leaving all of the J3 bonds 
satisfied. This structure has an additional 2-fold degen- 
eracy as columns can equally well be 'rows'. Again it is 
convenient to carry out a spin rotation on 'down' sites, 
yielding 
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with 



H = --{2J2 + J3) N + Ho + XV 



(5) 



Ho 
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(ly) 

<»i> 



QZ 



(2) 

^^E 

<ik> 



Sj St 
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^aE 

<il> 
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<ij> 
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Figure 4: Pliase diagram obtained from various series expan- 
sions. The dashed lines are first-order transitions, and the 
filled triangle is a triple point. 



Here the notations {Ix) and (ly) refer to nearest neighbor 
Ji bonds in the x and y direction in a plane, respectively. 
We have calculated the expansion to order A® inclusive, 
involving a total of 114650 distinct clusters with 4 types 
of bond. 



C. Results 

Series were computed and analysed for the cases J3 = 
0, 0.05, 0.1, 0.2, 0.5, 1.0 and a range of values of J2 in both 



the Neel and columnar phases. We set Ji = 1 through- 
out. 

Figure 2 shows estimates of the ground state energy 
per site versus J2 for two cases, J3 = (Fig. 2a) and 
J3 = 1.0 (Fig. 2b). The energy is evaluated by forming 
Fade approximants to the series and evaluating these at 
A = 1. Error bars, where shown, represent confidence 
limits, based on the degree of consistency between high- 
order approximants. 

As is apparent, the behaviour of the energy estimates is 
quite different in the two cases. For J3 = 0, i.e the J1 — J2 
square lattice model, the energy estimates in the inter- 
mediate range J2 = (0.4, 0.6) become erratic and the two 
curves, from the Neel and columnar series respectively, 
do not appear to join. However, for J3 = 1, the estimates 
are quite precise and the two curves cross near J2 ~ 0.55. 
The clear difference in slope of the two branches at the 
crossing point indicates a direct first-order transition be- 
tween the two phases. 

To analyse the magnetization series we first performed 
a Huse transformation (2l| to a new variable x ~ 1 — 
\/l — A to remove the square-root singularity at A = 1, 
expected from spin- wave theory. This procedure has been 
used in earlier work on the 2D antiferromagnet (21I . [2^ . 
Fade approximants to the new series were then evaluated 
at x = 1. 

The results are shown in Figure 3 and again display 
a striking difference between the two cases J3 = and 
J3 = 1. For J3 = we see clear evidence of an in- 
termediate non-magnetic phase, with Neel and columnar 
magnetizations vanishing near J2 ~ 0.4 and J2 ~ 0.6 re- 
spectively. The error bars are large near the transition 
points and it is not possible to say whether the transitions 
are first- or second-order. It is generally believed that 
the columnar-disordered transition is first-order while the 
Neel-disordered transition appears second-order but may, 
in fact, be weakly first-order [l^. On the other hand, for 
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Figure 5: Series estimates of the single- magnon dispersion relations in the Neel phase at J2/J1 = 0.0, J3/J1 = 0.5 (upper 
curve), and J2/J1 = 0.4, J3/J1 = 0.5 (lower curve). The dashed lines are linear spin-wave predictions. Values of q = k/vr along 
the path are shown at the bottom. 



J3 = 1 (Fig. 3b) it is clear that both magnetizations re- 
main finite throughout their phase. This is consistent 
with a direct first-order transition. 



Proceeding in this way with other values of J3, we 
construct a phase diagram as shown in Figure 4, which 
confirms the schematic phase diagram in Figure lb. As is 
apparent, the disordered phase narrows as the interplane 
coupling J3 is increased, and vanishes at the estimated 
point J2/J1 = 0.54 ± 0.03, J3/J1 = 0.16 ± 0.03. This 
value of J3 is a little lower than found previously [l6| . 



III. EXCITATIONS 

It is also of interest to investigate the spectrum of 
single-magnon excitations in the ordered phases and to 
see how these change on increasing the interlayer cou- 
pling J3. Standard series methods exist for calculating 
the excitation energy throughout the Brillouin zone [TtI , 
and we utilize these in the following. 



A. Neel Phase 

Series were obtained to order A^, based on 255196 dis- 
tinct clusters of up to 8 sites. In Fig. 5 we show disper- 
sion curves for two values of the frustrating interaction 
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Figure 6: Series estimates of the single-magnon dispersion relations in the columnar phase at (a) J2/J1 = 1.0, J3/J1 = 0.5, = 
(upper curve) and J2/J1 = 0.6, J3/J1 — 0.5, — (lower curve), and (b) J2/J1 — 1.0, J3/J1 = 0.5, fcz = 7r/2, (upper curve) 
and J2/J1 = 0.6, J3/J1 = 0.5,^2 = n/2 (lower curve). The dashes lines are linear spin-wave predictions. Values of q = k/yr 
along the path are shown at the bottom. 



J2 = 0.0,0.4 for J3 — 0.5, along various symmetry lines 
in the Brillouin zone, as shown. For comparison we also 
show the results of linear spin- wave theory. 
Several features are worth noting: 

(1) A very pronounced dip develops at wavevector 
(tt, 0,0) (and, by symmetry, also at (0, 7r,0)) with 
increasing J2. Linear spin- wave theory has the gap 
vanishing at J2 = 0.5 (see Appendix). 

(2) There is a distinct 'shoulder' at (0, 0, n/2), which is 
not seen in linear spin-wave theory. Indeed, along 



the whole line (0,0, fc^) the energy depends only 
weakly on J2. 

(3) Linear spin- wave theory reproduces the overall dis- 
persion curves reasonably well, but systematically 
underestimates the excitation energies. 

We note also that the series become very irregular near 
the Goldstone point k = 0, and the Fade approximants 
do not accurately show the vanishing of the gap at this 
point. This has been noted in previous work |23| . and an 
indirect way found to overcome this problem, which also 
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Figure 7: Series estimates of the energy gap e(k) at k = (0,7r,0) as a function of J2 at (a) — and (b) J3 = 0.5. Filled 
circles - Neel phase; filled triangles - columnar phase. The dashed lines are linear spin-wave predictions in each phase. 



allows calculation of the spin- wave velocity. However, we 
do not address this further in the present work. 



B. Columnar Phase 

Series were obtained to order A^, based on 490487 dis- 
tinct ckusters with up to 8 sites, and 4 possible bond 
types. In Figure[B]we show dispersion curves for J3 = 0.5 
and J2 = 1.0,0.6 along various symmetry liires in the 
Brillouin zone. The predictioirs of linear spiir-wave the- 
ory are shown as dashed lines. 

The most notable feature is the depressed excitation 
energy along the line (7r/2,7r,0) to (0,7r, 0). The disper- 
sion curve in this region flattens as J2 decreases from 
1.0 to 0.6, in both cases (a) and (b). According to liir- 
ear spin-wave theory this whole branch becomes soft at 
J2/J1 = 0.5, the classical transition point. Again liir- 
ear spin-wave theory tends to underestimate the excita- 
tion energy, but not as badly as for the Neel phase. The 
same difficulty with the Goldstone point k = 0, discussed 
above for the Neel phase, occurs here. 



C. The Energy Gap at (0,7r,0) 



of the gap is of some interest, and our results are shown 
in Figure 7, as a function of J2, for two values of the 
couphng ratio J3/J1 = 0, 0.5. 

Iir the case J3/J1 = 0.5 (Fig. [7)3) the gaps remain of 
substantial size right up to the direct first-order transi- 
tion, and the gaps appear to be of roughly equal magiri- 
tude in both phases at the transition. This is a further 
sign that the transition between the two ordered phases 
is first order at this value of J3. 

Figure [7^ shows the energy gap at J3 = 0. Note that 
whereas linear spin wave theory has the gap vanishing at 
all couplings in the columnar phase, it is in fact finite at 
larger J2, as given by a modified spin- wave theory and 
earlier series results [IjI- Note also that the energy gaps 
in both phases extrapolate to zero before reachirrg the 
classical transition point J2/J1 = 0.5, providing further 
evidence of a disordered intermediate phase when J3 = 0. 
Taken at face value, the results indicate that the gap re- 
mains finite at the estimated transition points J2 ~ 0.4 
and J2 — 0.6, which would indicate first order transi- 
tions. Naive extrapolations are not entirely reliable near 
a critical point, however, and a more precise check of this 
point would be of interest. A Dlog Fade analysis of the 
series in A gives no consistent estimates of the critical 
point iir this case. 



As the phase transition is approached, from either the 
Neel or columnar side, we expect increasiirg quantum 
fluctuations at wavevector k = (0, tt, 0), which should be 
manifest through a decreasing energy gap at this point. 
As mentioned above, linear spin-wave theory predicts a 
vanishing excitation energy at this k point at the clas- 
sical transition point J2/J1 = 0.5. Hence the variation 



IV. CONCLUSIONS 

We have used series expansion methods to study a sys- 
tem of stacked Ji — J2 layers with antiferromagnetic cou- 
pling J3 between the layers. This model is applicable 
to the layered materials Li2VOSi04 and Li2VOGe04, as 
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well as being of intrinsic theoretical interest. In agree- 
ment with earlier work of Schmalfuss et el. [l6j . we 
find that the disordered region of the phase diagram be- 
comes narrower as J3 is increased, and vanishes com- 
pletely at a triple point, beyond which there is a di- 
rect first-order transition between the Neel and columnar 
phases. We estimate the location of the triple-point as 
J2/J1 = 0.54 ± 0.03, J3/J1 = 0.16 ± 0.03, the J3 value 
being somewhat lower than in ref. [l6j . 

As well as mapping out the phase ground state phase 
diagram, we have computed magnon dispersion curves 
along symmetry lines in the Brillouin zone, in both the 
Neel and columnar phases, for various parameter values. 
This would allow a more critical evaluation of the va- 
lidity of the model for the materials mentioned above, 
when experimental inelastic neutron scattering results for 
magnon energies become available. 

After this work had been completed we became aware 
of two other recent studies of this model. Nunes et al. 
[2^ . using an effective field theory approach, concluded 
that the triple point (they refer to this as a critical end- 
point, assuming the Neel to disordered transition is sec- 
ond order) lies at J2/J1 = 0.56, J3/J1 — 0.67, values 
much higher than obtained in the present work and in 
other work [T6l [20j . Majumdar [2^ has carried out a 
spin- wave calculation to second order (l/S"^) and con- 
cludes that the intermediate phase remains present even 
at J3 = 1, which seems surprising in view of the other 
work referred to above. 

Finally we mention another recent study [l^ of a dif- 
ferent generalization of the 2-dimensional Ji — J2 model. 



with 

7fc = — (cos + cos ky), fii^ ~ cos k^ cos ky, Vf^ = cos k^ 

^ (A2) 
This is then diagonalized by a standard Bogoliubov 
transformation, yielding 

H = Eo + Y, M^Uk + BlBk) (A3) 

k 

with 

= 2S^Pl-Ql 

Pk = 2Ji + Ji - 2J2{1 - fik) 

Qk = 2Ji-ik + hvk (A4) 



in which the frustrating J2 interactions are included in 
all spatial directions. That model appears also to have 
an interesting phase diagram, and could be studied by 
series methods. 



V. ACKNOWLEDGEMENTS 

We are grateful for the computing resources provided 
by the Australian Partnership for Advanced Computing 
(APAC) National Facility. 



Appendix A: Linear Spin- Wave Theory 

In the Appendix we present the results of linear spin- 
wave theory (LSWT) for the model, for both the Neel and 
columnar phases. The results have been used in the main 
text to compare with the more accurate series results. 

The basic procedures of LSWT are well known so we 
will not include all details. 



1. Neel Phase 

There are two equivalent sublattices A, B and two sets 
of bosons describing spin deviations from the classical 
Neel state. After Fourier transformation the Hamiltonian 
is, to quadratic order 



(Al) 



The ground state energy is 
-Bo - -NS{S+l){2Ji - 2J2 + .h)+Y,^k (A5) 

k 

and the sublattice magnetization is 

We note that for small k 

Pk + Qk = 2i2Ji + J3)~^{J,+2J2){kl + kl) 

-ij3fc2 + ... 

Pk-Qk = {lji-J2){kl + kl) + y3kl + ---{A7) 



I 

H = -NS^{2Ji-2J2 + J3) + 2S{2Ji-2J2 + J3)Y,i4ak + blbk) 

k 

+ASJ1 lk{albl + akbk) + 4,3 J2 ^ /ifc(a^afc + &[.6fc) + 2SJ3 ^ MaWk + ^khk) 

I 
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so there is a linear Goldstone mode at k = 0, provided 
J2 < >^i/2. The phase becomes unstable for J2 > Ji/2. 

Another special case occurs at k = (tt, 0,0). Let k' = 
k — (tt, 0, 0) be small, then 

Pk+Qk = 2Ji-U2+'^k'l + Uk'l + k'l) + --- 
Pk-Qk = 2J1-4J2 + 2J3 + J2(fc''+fc'^) 

^|fc'^ + ... (A8) 



A Bogoliubov transformation then yields 

H = Eo + Y,^k{AlAk + BlBk) (AlO) 
fc 

with 

Eo = -NS{S + l){2J2 + .h)+^ujk (All) 

k 

and 

Pk = Ji cos ky + 2J2 + J3 

Qk = Ji COS kx + 2 J2 cos kx cos ky 

+J3Cosk^ (A12) 

The sublattice magnetization is again given by the for- 
mula ((A6| . 



[1] P. Chandra and B. Doucot, Phys. Rev. B 38, 9335 (1988). 
[2] E. Dagotto and A. Moreo, Phys. Rev. Lett. 63, 2148 
(1989). 

[3] H.J. Schulz and T.A.L. Ziman, Europhys. Lett. 18, 355 

(1992); H.J. Schulz, T.A.L. Ziman and D. Poilblanc, J. 

Phys. 16, 675 (1996). 
[4] J. Richter, Phys. Rev. B47, 5794 (1993); J. Richter, N.B. 

Ivanov and K. Retzlaff, Europhys. Lett. 25, 545 (1994). 
[5] J. Oitmaa and Zheng Weihong, Phys. Rev. B 54, 3022 

(1996). 

[6] R.F. Bishop, D.J.J. Farnell and J.B. Parkinson, Phys. 
Rev. B58, 6394 (1998). 

[7] R.R.P. Singh, Zheng Weihong, C.J. Hamer and J. Oit- 
maa, Phys. Rev. B60, 7278 (1999). 

[8] L. Capriotti and S. Sorella, Phys. Rev. Lett. 84, 3173 



Hence the gap at k' = vanishes at J2 = Ji/2, regardless 
of J3. corresponding to a phase transition. 



2. Columnar Phase 

There are again two sublattices, with the A and B sites 
forming successive columns. The Fourier transformed 
Hamiltonian is now 



(A9) 



For small k 



Pk-Qk = ^(Jl+2J2)fc2 + ^(2J2- Ji)fc2 

+ i Jsfc? + • • • (A13) 

which shows that again there is a Goldstone mode at k = 
0, which is stable for J2 > Ji/2, but becomes unstable 
at J2 = Ji/2. 

Another special case is k = (fc^,7r,0), when 



Pk - Qk = (2J2 - Ji)(l + cosfc,) + • • • . (A14) 

Hence uik vanishes along this line at J2 = Ji/2, corre- 
sponding again to the transition point. 



(2000). 

[9] L. Siurakshina, D. Ihle and R. Hayn, Phys. Rev. B64, 
104406 (2001). 

[10] O. P. Sushkov, J. Oitmaa and Z. Weihong, Phys. Rev. B 

63, 104420 (2001). 
[11] L. Capriotti, F. Becca, A. Parola and S. Sorella, Phys. 

Rev. Lett. 87, 097201 (2001). 
[12] L. Capriotti, F. Becca, A., Parola, A., and S. Sorella, 

Phys. Rev. B 67, 212402 (2003). 
[13] R.R.P. Singh, Weihong Zheng, J. Oitmaa, O.P. Sushkov 

and C.J. Hamer, Phys. Rev. Lett. 91, 017201 (2003). 
[14] T. Roscilde, A. Feiguin, A.L. Chernyshev, S. Liu and S. 

Haas, Phys. Rev. Lett. 93, 017203 (2004). 
[15] J. Sirker, Z. Weihong, O. P. Sushkov and J. Oitmaa, 

Phys. Rev. B 73, 184420 (2006). 



J 



H = -NS^{2J2 + J3) + 2SY,iJicosky + 2J2 + J3){alak + blbk) 

k 

+2S ^(Ji cosfcj: 4- 2J2/ifc + J3'^k){albl + aubk) 

k 

I 



10 



[16] D. Schmalfuss, R. Darradi, J. Richter, J. Schulenberg and 
D. Ihle, Phys. Rev. Lett. 97, 157201 (2006). 

[17] J. Oitmaa, C. J. Hamer and Z. Weihong, Series Ex- 
pansion Methods for Strongly Interacting Lattice Models, 
Cambridge, (2006) 

[18[ H. Rosner, R. R. P. Singh, W. H. Zheng, J. Oitmaa, 
S.-L. Drechsler and W. E. Pickett, Phys. Rev. Lett. 88, 
186405 (2002); H. Rosner, R. R. P. Singh, W. H. Zheng, 
J. Oitmaa, and W. E. Pickett, Phys. Rev. B 67, 014416 
(2003). 

[19[ G. S. Uhrig, M. Holt, J. Oitmaa, O. P. Sushkov, and R. 

R. P. Singh, Phys. Rev. B 79, 092416 (2009). 
[20[ M. Holt, O.P. Sushkov, D. Stanek and G.S. Uhrig, 



preprint [i^iv:1010.555lV l (2010) 
[21[ D. A. Huse, Phys. Rev. B 37, 2380 (1988). 
[22[ Zheng Weihong, J. Oitmaa and C.J. Hamer, Phys. Rev. 

B 43, 8321 (1991). 
[23[ R.R.P. Singh, Phys. Rev. B47, 12337 (1993); R.R.P. 

Singh and M.P. Gelfand, ihid. 52, R15695 (1995). 
[24[ W. A. Nunes, J. Ricardo de Sousa, J. Roberto Viana and 

J. Richter, J. Phys. Cond. Mat. 22, 146004 (2010). 
[25[ K. Majumdar, J. Phys. Cond. Mat. 23, 046001 (2011). 
[26| K. Majumdar and T. Datta, J. Stat. Phys. 139, 714 

(2010). 



